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Abstract: We consider QCD factorization between hard and soft subprocesses in inclusive 
reactions where the momentum fraction x of one parton approaches unity as the hard scale 
Q 2 — ► 00, such that Q 2 (l — x) is fixed. In this "BB limit" the entire (multi-parton) 
Fock state containing the high x parton is coherent with the hard subprocess. The soft 
contribution is given by a forward multiparton matrix element. The BB limit corresponds to 
a fixed (large or small) missing mass and is thus closely connected to exclusive production. 
We analyze the Drell-Yan process h + N — > 7* + X in detail, explaining why the virtual 
photon is longitudinally polarized for h = ir and transversely polarized for h = p. The BB 
limit may be relevant also for other phenomena observed at high x, such as the large single 
spin asymmetries of pp — > X and in pp^ — > irX. 
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1. The BB limit 



Data on hard inclusive processes has been successfully analyzed assuming QCD factoriza- 
tion between a hard subprocess and universal soft matrix elements (parton distributions). 
Formally one considers the leading contributions in the Bj limit where a hard scale Q 2 — > oo 
while the momentum fractions Xi carried by the active partons (one in each hadron) are 
held fixed. The higher twist corrections are power suppressed in the hard scale but gener- 
ically increase as Xi — ► 1. Thus the effective expansion parameter is 1 /[Q 2 (l — (see [1] 
for a recent phenomenological analysis of eN — ► eX (DIS)). 

Data at high x has features which differ qualitatively from the leading twist contri- 
bution. A striking example is the polarization of the virtual photon in ttN — ► fj, + fj,~X 
which changes from transverse to longitudinal when the photon carries a momentum frac- 
tion xp > 0.6 of the pion beam [2]. This requires that at least one of the annihilating 
quarks in the subprocess qq — > 7*(Q 2 ) goes off-shell by an amount commensurate with 
the photon virtuality Q 2 . As pointed out by Berger and Brodsky [3] the longitudinally 
polarized photon is coherent with both valence quarks in the pion. Hence the dominant 
contribution is of higher twist even though the subprocess is hard. An analogous change 
of J /ip polarization at high xp was observed in irN — > J /tp + X [4]. The inapplicability of 
the twist expansion in QCD at high x was discussed in [5] . 

In this paper we consider QCD factorization in a limit where x — > 1 as Q 2 — > oo. Here 
x can denote either the momentum fraction Xj of a fast quark i or equivalently the xp of a 
final state particle. The Drell-Yan data suggests that the helicity of the pion is transferred 
to the virtual photon at high x, implying that the photon is coherent with an entire Fock 
state of the pion. The life-time of a Fock state is inversely proportional to AE, the energy 
difference between the pion and the Fock state. At high pion momentum p, 

2pAE^ml-^2 p2i± + m2 > (1.1) 

Xi 

i 

The x ~ 1 quark which annihilates into the virtual photon contributes only little to the 
energy difference (1.1). The life-time of the Fock state is determined by the stopped partons 
with Xi oc 1 — x. If 1 — x ~ Aq CD /Q 2 we get 2pAE ~ Q 2 (we take pi± ~ Kqqd, a generic 
soft QCD scale). Such Fock states have life-times similar to that of the virtual photon, 
ensuring their coherence with the hard process. This motivates us to consider the 

BB limit : Q 2 -> oo at fixed Q 2 {1 - x) (1.2) 

which recognizes [6] the early observations of Berger and Brodsky [3] . 

The understanding of QCD factorization in the BB limit is facilitated by an analogy 
with ordinary DIS {ep — > eX). In the rest frame of the target DIS may be viewed [7] 
as proceeding through a splitting of the virtual photon into a quark pair, 7*(Q 2 ) — ► qq. 
The quark carries nearly all of the photon energy v and forms the current jet, whereas the 
antiquark carries a finite momentum in the target rest frame, even as v = Q 2 /2m p XB — ► oo. 
The qq Fock state of the photon is analogous to the pion Fock state we just discussed: A 
fast quark with momentum fraction x q ~ 1 and an antiquark with Xq ~ Aq CD /Q 2 . The 
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fact that the momentum Xqf of the antiquark remains finite in the Bj limit allows the 
alternative (standard) interpretation of the antiquark as a quark constituent of the proton. 
The DIS cross section can then be equivalently ascribed to the scattering of the qq Fock 
state on the target or to the probability of finding a quark in the target wave function. The 
latter interpretation corresponds to the parton distributions obtained in QCD factorization. 

In the following we shall analyze factorization of the Drell-Yan process in the BB limit, 
interpreting stopped partons in Fock states of the beam as (anti-)parton constituents of 
the target. We do not consider rescattering effects, i.e., gauge links between the quark 
fields [8]. 



2. Kinematics of ir + N -> 7* + X at fixed (1 - x F )Q 2 

Our notation for the ir + N — > 7* + X kinematics is indicated in the Feynman diagram of 
Fig. 1. The pion, nucleon and photon momenta are in Light-Front (LF) notation (k = 
(k + ,k~,k ± ) with fc± = k°±k 3 ) 



Pion : 
Nucleon : 
Photon : 



(0,k-,0±) 

(p + ,m 2 N /p A 

(Q 2 Ar,T 



,o±) 

q±) (2.1) 



We neglect the mass of the pion as we take 
its momentum k~ — ► 00. We work in the 
target rest frame, p + = p~ = mjv, thus 
s = (k + p) 2 ~ rriNk~ — ► 00. In the BB 



limit xf = 
ing also 

x B 



q jk — > 1 at fixed q ± , keep- 



q 



±- = 
p+ 

XMP~ 



Q 2 

s 



fixed 



Q 2 {l-x F ) 



(2.2) 
fixed 




GPD 



Figure 1: Feynman diagram contributing to the 
Td^N — > 7*^0 amplitude. The arrows indicate 
momentum directions; the ?/, are space-time po- 
sitions of the vertices and A, B denote color in- 
dices. The pion contributes via its ud distribution 
amplitude 4> n and the GPD denotes a transition 
(N — > X) Generalized Parton Distribution. 



Since the four-momentum transfer k — q to the target system is fixed the invariant mass 
Mx of the hadronic state X is finite in the BB limit, 



M 



x 



(k + p - q) 2 ~ (1 - xb)[s(1 - xf) + m 2 N ] 



(2.3) 



and may be small or large depending on the value of s(l — xp) = Q 2 (l — xf)/xb- For 
X = N we have the exclusive Drell-Yan process ttN — ► 7*^ [10] as well as the time-reversed 
version of 7*iV — ► tt + N (Deeply Exclusive Meson Production) which are well-known to be 
described by Generalized Parton Distributions [11] as indicated in Fig. 1. Factorization 
[12] applies equally to the transition GPD's (X ^ N) [13, 14] since Mx *C Q. In section 3 
we merely sketch the derivation of the factorized amplitude. 

In section 4 we consider the factorization of the inclusive Drell-Yan cross section 
a(n + N — > /j,[j,X) through a completeness sum over the states X. As in the standard leading 
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twist limit of the Drell-Yan process, this includes integrating over the (small) transverse 
momentum q± of the virtual photon, which arises from the intrinsic momentum of the 
partons in the GPD. High transverse momenta are suppressed at fixed Mx since the GPD 
is expected to fall off oc similarly to nucleon transition form factors. The integration 

over can be done at fixed Mx by varying (k — q)~ without affecting the hard subprocess 
at leading order. 

We anticipate that the pion will contribute through its distribution amplitude (^^(z) 
and parametrize the momenta of its valence quarks as 

ki = (0, zk~ , 

k 2 = (0,(1- z)k~,-k ± ) (2.4) 

We take the transverse momenta to be limited, k\ <C Q 2 ■ Thus we neglect the perturbative 
tail of the wave function which arises from gluon exchange and gives rise to the logarithmic 
Q 2 evolution of the distribution amplitude [9] . kf and k^ vanish in the high energy limit 
and do not contribute at leading order. 

In the BB limit (1.2) one of the quarks in the pion (the u-quark in Fig. 1) transfers 
nearly all its momentum to the other quark, such that the photon carries nearly all of the 
pion momentum. The stopped u-quark is left with a finite momentum in the target rest 
frame and should be connected to the target wave function as discussed above for DIS. 
Consequently, as k~ — > oo, 

x = if /p fixed 
xb + x = I2 /p + fixed (2-5) 

A large longitudinal momentum is transferred through the gluon (q\) and quark (q 2 ) prop- 
agators whose virtualities are of O (Q 2 ) , 



q\ ~ — zk I 
q 2 c -k-lf 



cxQ 2 = x B s (2.6) 



This implies that the hard interactions occur at nearly the same LF times (yf, y£, yt = 
O (l/k~)) and vanishing transverse separations (yu_, y 2 ±, 2/3_i_ = O (1/Q)). 

Since q 2 and q\ are independent of l~ and li± the soft target matrix element will 
be integrated over these components (with the kinematic constraint l\ — l 2 = k — q). In 
evaluating the hard, perturbative subprocess it suffices to take 

h = (xp+,0-,O ± ) 

h = ((x + x B )p + , 0~,0 ± ) (2.7) 

and analogously only the minus components in (2.4) contribute. Thus both quarks effec- 
tively reverse their direction of motion along the z-axis. 
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3. Drell-Yan amplitude in the BB limit 



The general expression for the scat- 
tering amplitude is 

iT(-K + N -» j* x X){2ir) 4 5 4 {k +p-q-p) 



( 7 ^|T|exp[- 



dtH T (t)] }\ttN) 



(3.1) 




The expression for a diagram such as in 
Fig. 1 is obtained by expanding the hard 
vertices of the interaction Hamiltonian Hj, 
connecting them by perturbative propaga- 
tors and retaining the leading contribution 
in the BB limit (2.2). Equivalently, we 

may simply note that due to the large momenta of 0{k~) flowing from y\ 
the corresponding LF time differences yf — yf and yf — yf are of 0(l/k~ 



Figure 2: Feynman diagram contributing to 
T u (tt + N -> j*X). Notations as in Fig. 1. 



2/2 -> V3 

Similarly, 

the squared transverse separations vanish as the inverse virtualities (2.6): \y2± — Ui±\ and 
|y 3 _i_ — 3/ 2 _lI are °f O 0-/Q)- The diagram may then be evaluated using standard Feynman 
rules with the prescription (A. 14) (see Appendix A) for the pion. The quarks l\ and I2 are 
treated as external particles and their free wave functions are replaced according to 



u(h) 
u{h) 



2vr 2 



(X(p')\il) u (y 1 )eyLY>{-\iy 1 if 



(3.2) 



Our conventions for the wave functions and operators are explained in detail in Appendix A. 
Adding the diagram where the gluon vertex yi in Fig. 1 is on the I2 line we find for the 
amplitude where a photon of helicity A = is emitted from the <i-quark, 



r d (7r+iV - 7 £X) = e d - 



-ieg 2 Cf 



(3.3) 



/ dzdl tMz) I d y-e-^ l V\x{p')\Uyi)i + iMmip)) yt=m 



z2tt(IJ 



le) 



where 4> n (z) is the pion distribution amplitude. The amplitude for transversely polarized 
photons is suppressed by a factor 1/Q. We give an intuitive explanation of this in section 5. 

The contribution from the diagrams where the photon is emitted from the u-quark as 
in Fig. 2 is obtained in a similar way. Defining 



C(x B ,x) 



dz (ft-xtz) 



1 — zxs + x + ie 



+ 



ed 



1 



z x 



IE 



(3.4) 
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the Drell-Yan amplitude for longitudinal photons is 

,2 



T^N ,IX) =T U + T d = / dx C(x B , X ) 

(3.5) 

x I dyre-^^ /2 (X(p')|^(yi)7 + 75V'd(0)|iV(p)) 



yt=yix=° 

where x and xb are defined in (2.5). For X = N we may verify that the amplitude indeed 
corresponds to the usual expression for deeply exclusive pion production [12]. For general 
states X the matrix element is a "transition" GPD [13, 14]. 
The full amplitude including the muon vertex is 

T(ttN n + n~X) = T(itN rtX)-^u( qi ,s)(-ie)t (q)v(q 2 , -s>) 

= T(kN -> 1 lX) ie ^5 ss , (3.6) 

where is the polar angle of the muon momentum in the muon pair rest frame and s, s' 
are the muon helicities. 



4. The Inclusive Drell-Yan cross section 

The inclusive cross section is obtained by summing over all final states X, 

o{k+N^ 1 IX) = ± s Y, j ^0^\T(tt + N ^rtX)\\2n)W(k + P - q - P ') (4.1) 
where 1S defined by the completeness relation, 

£ \X)(X\ ee £ / J] J 2 J^\P 1 , • • ■ , Pn )( Pl , ..., Pn \ = l (4.2) 

X n=0 J i=l V 7 * 

and p' = X^jPi i s the total momentum of the state X. 

The momentum conserving (^-functions in (4.1) constrain p' and thus allow only a 
subset of the states X in the completeness sum (4.2). The restriction on p\ is avoided by 
summing over all transverse momenta q± of the virtual photon. As mentioned in section 
2, large values of q± do not contribute significantly to the sum due to the suppression 
provided by the GPD's. On the other hand, p'+ and p'~ are fixed by a measurement of 
xb and xf- Hence we need to incorporate the corresponding 5- functions in the (hermitian 
conjugate) matrix element: 

(N(p) 1^,(0)7+76 My2)\X(p'))2(27T) 2 5(p + -q + -p' + )5(k~ +p- - q~ - p'~) 

= \ J dy+ dy^ (iV(p)|V^(y 3 )7 + 75 MV2 + V3)\Xtf)) exp [iy 3 ■ (k - q)] (4.3) 

where y 3 j_ = and we used translation invariance in the matrix element. The full complete- 
ness sum is allowed using the matrix element in (4.3), since states X which do not conserve 
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Figure 3: Pictorial representation of the forward multiparton distribution fdu/p( x B, %M] x, x') 
given in Eq. (4.4). 



momentum will not contribute after the 2/3-integrations. The inclusive cross section (4.1) 
will thus be given by the multiparton distribution shown in Fig. 3, 



fdu/ P {xB,x M ;x,x') = (4.4) 

x(N(p)\^d(y3h + l5^u(y2 + 2/3) V^G/l)7 + 75 ^d(0)\N(p)) y . ±=0 . y+=y+=0 

where xb = Q + /p + , x = if /p + , x' = l'i /p + and the scaled ' — ' momentum transferred to 
the inclusive system is denoted xm = k~(l — xp)/p~ ■ The inclusive mass Mx is given by 

xm as 

Ml = 77^(1 - x fl )(l + x M ) - q± (4.5) 
The kinematic range of xm hi the BB limit is thus 

x B + q\/m 2 N 

: < x M < 00 (4.6) 

1 - XB 

The shift by introduced in (4.3) between the fields in the matrix elements of T 
and in is conjugate to the momentum transfer k — q between the hard vertex and the 
target. The MPD (4.4) differs from the higher twist distributions discussed by Jaffe [15] 
through its dependence on the finite LF time difference between the fields. Also the 
standard leading twist PDF's are evaluated at y% = since the inclusive system X carries 
an asymptotically large '— ' momentum in the Bjorken limit. 

In the BB limit the the momentum transfer xmP~ = k~(l — xp) = {l\ — Z2) - ~ I/2/3 
to the target is kept finite. We may, however, consider the case where this transfer is 
nevertheless large compared to p~ . Since the quark lines l\ and I2 connect to the non- 
perturbative matrix element their virtualities if, 1% should remain limited, thus if ~ 
Aq CD /1^ and similarly for if . Given that if — if = q + = xbP + is fixed, either l^ or —l^ 
may be large, while the other has to be of O (p~). If V[ is large then if ~ and if ~ q + . 
In this kinematics the u-quark will hadronize independently into a final state jet [3] and the 
target MPD (4.4) should reduce to a d-quark PDF. We demonstrate this in Appendix B. 
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Using the expression (3.5) of the scattering amplitude in the cross section (4.1) and 
incorporating the longitudinal (^-function constraints in T' as in (4.3) the Drell-Yan cross 
section in the BB limit can be expressed in terms of the double parton distribution (4.4) 
as 

^/t2 1lX ^ = CF \at f dxdx' C(x B ,x)C*(x B ,x') f du / p (x B ,x M ;x,x') 

uJVL-£ S [1 — Xb)^c J 

(4.7) 

Including the muon pair observables as in (3.6) the differential cross section becomes 

da(ir + N -> ppX) _ 2 sin 2 6 da(ir + N -► ~f* L X) 
dx B dtt^ dM\ ~ (4vr) 3 Xij ~dM\ ^ 



5. Photon helicity for pion and proton induced Drell-Yan processes 

The fact that the photon is coherent 
with an entire pion Fock state makes it 
natural that it also carries the helicity (A = 
0) of the pion. On the other hand, this ar- 
gument does not suffice to determine the 
photon polarization in pN — > j*X, since 
there is a minimal helicity flip | AA| = ^ f° r 
both transverse and longitudinal photons. 
The photon was found to be transversely 
polarized in a calculation of the exclusive 
pN — > 7*7r process [14]. As we shall see, 
the same result is obtained for the proton 
induced inclusive Drell-Yan process in the 
BB limit. 

The helicity systematics follows in a straightforward way from three facts: 

(i) The hard interactions conserve quark helicity up to corrections of O (m q /Q); 

(ii) Since all transverse momenta q± are limited, orbital contributions L z ~ O (q±/Q); 

(iii) Angular momentum J z ~ S z is conserved. 




Figure 4: The double red arrows indicate the spin 
directions of the particles. All momenta are in the 
indirection as shown in the parentheses. The it- 
quark propagates into the GPD (not shown) while 
the d-quark propagates out of it. 



The helicities are then obtained by simple addition. In Fig. 4 the J z components 
are indicated by double arrows in the case x > where the u-quark propagates (with 
positive energy) into the GPD while the d-quark propagates out of it. They follow from 
helicity conservation for the quark lines, taking into account the direction of longitudinal 
momentum of the particles (indicated by ± for motion in the ±z direction). E.g., the 
dominant spin component S z = — 1 of the gluon may be verified directly by expressing its 
propagator as a sum over helicities as in (A. 3). The quark propagators can analogously be 
expanded in terms of the spinors (A. 10). 
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The photon polarization in the case of a J = ^ beam particle, such as a nucleon, may 
be similarly deduced. The analysis is especially simple if we model the nucleon as scalar 
diquark-quark bound state. Setting S z (u) = in Fig. 4 gives S z (g) = and consequently 
S z (l*) = 1. 

We have calculated the relevant diagrams (Fig. 5) also for three-quark Fock states of 
a nucleon (see [14] for a corresponding calculation of the exclusive pN — ► 7*7r process). 
Ignoring flavor, the hard subprocess involves nine topologically different diagrams at tree 
level. See Fig. 6 for momentum and helicity labels. We parametrize the momenta in 
analogy to (2.4), (2.7) as 



h = 


(0+,zi£T,fc 1± ) 


k 2 = 


(0 + ,z 2 k-,k 2± ) 


k 3 = 


(0 + , (1 - z\ - z 2 )k~,-k l± - k 2 ±) 


li = 


{x lP + ,0-,£ 1± ) 


h = 


{x 2P + ,0-,£ 2± ) 


is = 


{-{x B + xi + x 2 )p + , 0", -£i± - £ 2 ± 


1 = 


{x B p + ,x F k~ ,0_l) 



where all transverse momenta (~ Aqcd) are negligible at leading order in 1/Q when 
evaluating the hard subprocess. 

We illustrate the method by calculating the diagram of Fig. 6 in the BB limit in 




7. 8. 9. 



Figure 5: Nine topologically different subprocess diagrams in pp — * 7*(xf) + X at large xf- 
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Feynman gauge. We suppress all color indices and use the light front spinors of (A. 10) for 
the proton valence quarks. Helicity conservation along the quark lines fixes £1,2,3 — s i,2,3> 
leaving si,s 2 , S3, and A as parameters. We find that four helicity amplitudes M^ 



receive contributions at leading order (M* 



i/Q 



3\. 



4\/2e q g i ^/x 2 (x 1 + x 2 + x B )z 2 (l 



z\ - z 2 



(p+k ) 3 / 2 y/xiZi(xi + x 2 ) 2 (zi + z 2 ) 



~ Ml 



~ Ml 



(5.2) 



where the lower indices ± stand for ±^ and e q is the quark charge at the photon vertex. The 
virtual photon is transverse in the leading amplitudes. This result can also be deduced by 
using helicity and angular momentum conservation as explained above. The z-components 
of spins are marked with red double arrows in Fig. 6 for the amplitude .M^ 1 . + . Notice 
that the spin directions match with those of the scalar diquark model above. 

We have checked by explicit calculation 
that the above reasoning applies similarly to 
all the diagrams of Fig. 5. We find that the fol- 
lowing helicity combinations are leading (fixing 
A = +l) : 



Diagram: 
1,4,8,9 
5 
6 
7 



Helicities: 
(si,s 2 ,s 3 ) 
(si,s 2 ,s 3 ) 
(si,s 2 ,s 3 ) 
(si,s 2 ,s 3 ) 



(+,-+), 

(-+,+) 
(+,+,-) 
(+,-+) 



(-+,+) 



(5.3) 




Figure 6: Definitions of momenta (fcj, £i) and 
helicities (sj,ij) in diagram 1 of Fig. 5. The 
red arrows indicate S z in units of when 
Sj = — S2 = S3 = \ and the photon has helic- 
ity A = 1. The second gluon carries S z = 0. 



Diagrams 2 and 3 do not contribute in Feyn- 
man gauge. Contributions with A = — 1 are 
obtained from (5.3) by flipping all quark helic- 
ities, and amplitudes with A = are subleading 
as expected. 

Interestingly, the (+, +, — ) amplitude gets 
a contribution only from diagram 6 of Fig. 5. 

This can be understood in terms of helicity conservation since the rules (i)-(iii) imply that 
for all other diagrams of Fig. 5 except diagram 6 either a gluon has \S Z \ = 2 or a quark 
has \S Z \ = 3/2. Note that diagram 6 is gauge invariant at leading order in 1/Q since both 
gluons are transverse. 

It is also straightforward to check that amplitudes where a pion (proton) projectile 
produces a transverse (longitudinal) virtual photon are suppressed by orbital factors of 
\£i,2±\/Q or \ki )2 ±\/Q as expected from rule (ii) above. In the suppressed amplitudes the 
projectile wave function is strongly weighted at the endpoints z = 0, 1. E.g., the pion 
distribution amplitude is weighted by 1/z 2 for z — > in the transverse photon amplitude, 
instead of by 1/z as in (3.4). This endangers the convergence of the z-integral and is related 
to the difficulty of proving factorization for 7^, N — > irN [12] . 
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6. Discussion 



The remarkable equality of the inclusive 7* TV" — > X DIS data at high Q 2 with the exclusive 
resonance contributions at lower Q 2 [16] points to an intriguing relation between inclusive 
and exclusive processes. This Bloom-Gilman duality indicates that two mathematical limits 
simultaneously describe the data: Q 2 — ► 00 at fixed xp and at fixed hadronic (resonance) 
mass W. The latter limit is equivalent to the BB limit of keeping Q 2 (l — xb) fixed, some 
aspects of which we have explored in this paper. 

We focussed on the Drell-Yan process ttN — > 7*X at fixed Q 2 (l — xp) for the virtual 
photon. The xp — > 1 limit forces a single quark to carry nearly all the pion momentum. 
Knowing the "end-point" behaviour of the pion wave function is important in the analysis of 
its exclusive form factor. The DY data [2] shows that the photon is longitudinally polarized 
at high xf, implying that the annihilating quark(s) are far off-shell. This dynamics was 
studied by Berger and Brodsky [3], who showed that the hard subprocess is different from 
the standard qq — ► 7* one and, as in exclusive processes, is coherent with the entire pion 
Fock state. 

The theoretical and experimental evidence of the dominance of a hard subprocess 
which is suppressed (i.e., of higher twist) in the usual limit (Q 2 — ► 00 at fixed xf) raised 
the question whether one can define another limit where this subprocess would be a leading 
contribution and could be factorized from the soft dynamics. As we discussed in section 1, 
the requirement of coherence naturally selects the BB limit of fixed Q 2 (l — xp). Since the 
mass of the inclusive system X is fixed in this limit the existing factorization proofs for 
deeply exclusive meson production 7*iV — ► ttN [12] apply. The derivation is valid for any 
fixed mass state X in irN — > 7*X, which builds the inclusive Drell-Yan process. The elastic 
process irN — > 7*iV was already considered in [10]. For I/iV the soft dynamics involves 
"Transition Generalized Parton Distributions", which have also been studied previously 
[13, 14]. 

We found that irN — > 7* A" differs essentially from pN — > 7* A in the BB limit since 
the photon is longitudinally polarized in irN and transversally in pN. Hence the transition 
from the Bj to the BB limit can be directly observed in ttN through the change of photon 
polarization with increasing xp. In pN there is (as expected) no evidence for a change 
of polarization from transverse to longitudinal, and the two limits may coexist in a range 
of xf as suggested by duality [16]. This may be connected to the observation that the 
azimuthal distribution of the muon pair in pN — ► /i + /i~A is consistent with the standard 
QCD analysis [17] while this is not the case for ttN — > fj> + /i~X [18]. 

The BB limit seems appropriate for understanding the large single spin asymmetries 
(SSA) observed in hadron scattering at high xp and large transverse momenta [6]. The 
SSA requires quark helicity flip and a dynamical phase, both of which are suppressed in 
hard subprocesses. The coherence between partons of high and low x allows the helicity 
flip and phase to occur in the soft part of the amplitude where they are not suppressed. 

The multiparton distribution (4.4) that describes the soft dynamics of the Drell-Yan 
process in the BB limit is a forward target matrix element with four quark fields, which 
would be of higher twist in the usual Bj limit of fixed xp. It differs from the MPD's 
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studied by Jaffe [15] in that the two pairs of quark fields originating from the amplitude 
and its complex conjugate are evaluated at a finite LF time difference y^~. This is a 
consequence of the finite (in the target rest frame) 'minus' momentum (1 + xm)p~ of the 
inclusive system and means that the ordering of the quark fields in the matrix element is 
significant. For xm — > oo we demonstrated that the main contribution to the MPD is from 
a contraction of a pair of quark fields, which reduces the MPD to a standard leading twist 
parton distribution. 
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Appendix 

A. The virtual photon and pion wave functions 

We use LF polarization vectors [19] for the virtual photon with helicity A, 

£x(q) = e\~ e ^ L n (A.l) 



with 

V2 



e±i = — ~ n = (0,0,±l,i) n= (0,2,0,0) 



(A.2) 



eo(q) = —^= h= (2,0,0,0) 



which satisfy 



£ el{k)e\{kY = -g^+ K -^- (A.3) 



k 2 
A=±1,0 



The pion valence state is defined by its wave function 'I'fc, 

br(fc)) = y / ^i^( fei) b i {kl)d t (fc2) | 0) (A . 4) 

1 V " V2Wc J V 4 l fc ill fc 2|(2vr) 3 S ' M ' ~ sAK n ' V ' 

where &2 = k — k\. The sum is over the helicity s = =t^ and color A (N c = 3). With 

{b s , A (kx), ^^(fe'i)} = (2vr) 3 2|fc 1 |5 3 (fc 1 - tiJSnfSAA, (A.5) 
this gives the normalization 

(vr(fc')lvr(fc)) = (2vr) 3 2|fc|5 3 (fc - k>) J ^*^ \* k ( hl )\ a (A.6) 
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Due to the large virtuality of the gluon propagator q\ in Fig. 1 only transversally compact 
valence Fock states contribute at leading order, involving the pion distribution amplitude 

Mz) = I ^y^ k[kl) (A - 7) 

where k± = zk~ as in (2.4). With f n ~ 93 MeV the normalization is 

1 f 

dz4>*{z) = %= (A.8) 



2JN r . 



o 



Since we took the pion to move along the — z direction in (2.1) we define the LF spinors 
[19] in the quark operator 

Mv) = f 2 \p\(Lr £ [ b ^ u °^ e ~ ip ' y + 4(p)Mpy p - y ] (a.9) 

to be well-defined for k + = 0, 

VP 

v(p,s) = —=fx(-s) (A.10) 

Vp 



where X {\) = (1,0,1,0)^/^2 and x (-i) = (0, 1, 0, -1) T /V2. 

At lowest order the Bethe-Salpeter wave function of the pion is 

*f (2/1,2/2) = (0\^(y 2 )^(yi)\TT(k)} (A.ll) 

-s)u a (k 1 ,s)exp(-iki ■ yi - ik 2 ■ y 2 ) 



V-/ n x,_i f d 3 fci^fc(fci) _ „ 



V2iV c 

For the LF spinors defined in (A. 10), 



^(-l)-5^(fe 2 , -^ttaCfci, 8) = -jL= ^(-l) S -5(^) Qa/X a'(s)(7 x)^( S )(fe)^ 



A; 1 fe 2 



[^^75^2] a/3 - y\fh k 2 [#75] a/3 (A.12) 



2^ k 1 k 2 

where the last equality applies in the limit — > 00, using ^ ~ Thus 

^--hx, (2/1, 2/2) = 2 ^^ [^75] Q /jfc" y ^fe(fcf = -|2fe _ ,fe_ L i)exp(—jA;ryi-ifc2-2/2) 

(A.13) 

With kf = k% = as in (2.4) the BS wave function ^-^^ (yi ,2/2) is independent of y\ , y 2 ■ 
Thus the wave function is the same whether evaluated at yj~ = y 2 (the LF wave function 
for a particle moving in the — z direction) or at equal ordinary time, y\—y 2 = ut ~ V\ • 
This expresses the well-known fact that equal time wave functions are the same as LF wave 
functions in the infinite momentum frame. 
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In the BB limit we have y%\_ = (i = 1, 2), hence the pion enters through its distribution 
amplitude (A. 7). In analogy to (3.2) we can write an effective Feynman rule for the pion 

up(h)v a (k 2 ) - J exp (|tfcTy+) (0|^|(0)^(y 1 )|7r(fc)) 2/r=2/ix=0 (A.14) 

which will be multiplied by = z/c~-dependent propagators, cf. (2.6). 



B. The limit of large xuV — k (1 ~~ x f] 



In the BB limit of the DY process ir + N — > "y*X the virtual photon carries nearly all 
{xf — ► 1) of the pion momentum, leaving a finite (in the target rest frame) transfer 
xmP~ = k~(l — xf) to the target system. In this Appendix we consider the case xm S> 1, 
when the momentum transfer (Zi — l2)~ is large compared to p~ '. Since the quark lines 
l\, I2 attach to the non-perturbative GPD their virtualities I 2 = — l 2 j_ should remain of 
O {^qcd) • Given that (I2 — h) + = xbP + is fixed either or — l 2 , but not both, can grow 
large. When Zj~ is large the w(/i)-quark hadronizes independently of the target remnants, 
similarly to the struck quark in the standard Bjorken limit of the DY process. Here we 
show that the target MPD in the expression (4.7) of the DY cross section reduces to the 



ci-quark PDF as xmP 



00, and to the u-quark PDF as xmP 



00. 



(i) fdu/p -> fd/p as xm -> 00 with x M p ^ h 

According to the expression (4.4) of the MPD the LF time difference between the quark 
fields y^ ~ 1/xmP~ — ► as xm — > 00. The dynamics then becomes light-cone dominated 
and contractions of the u- or <i-quark fields in the MPD dominate. A contraction turns 
the MPD into a standard PDF in which the field ordering is irrelevant, i.e., only one 
contraction is possible. 

The contribution where the «-quark has large momentum and is treated as a free final 
state particle (l± — > 00, IJ ~ ^-qcd/^i) corresponds to contracting the u-fields in fdu/p of 
(4.4). Using 



(0|^(x)^(0)|0> 



dH 



-il-x 



(2vr) 3 2|Z| 



i°=\i\ 



dl + d 2 l_ 
(2tt) 3 2Z- 



-e{i + )e- ilx t 



l-=l[ 1/1+ 



(B.l) 



we have 



4(4vr) 3 
dl+ d 2 l 



x exp 



(2tt) 3 2/+ 



6(1+) 2l + (N(p)\My 3 ) 7 + M0)\N(p)} 
I 2 



2 il + (y 2 +y 3 -y x )- 5*^1/3 



(B.2) 
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Integrating over the transverse momenta, 

2ttU + 



d 2 l± exp 



vt ~ i£ 



(B.3) 



Given the small effective range of yt ~ (1 — x^)) we may neglect the y^-dependence 

of the matrix element in (B.2) (which is regular on the light-cone). We get 

, exr>\},iytk^(l — xp)} 

dyt 2 3 , — = 2vri (B.4) 

yj-ie 

The integrals over and y^ give 6- functions constraining if = l[ + and l + = if. Noting 
that if = q + + if we have 

f d u/ P = S(lf ~ ^ 0(lt)fd/ P (lf/p + ) (B.5) 
where the d-quark PDF is 

f d / P (lf/p + ) = dy^ exp(-£n,3 1+) (N(p)\My3h + M0)\N(p)) (B.6) 

The dominant contribution to the cross section (4.7) comes from the region of low 
if = xp + due to \C(xb,x)\ 2 oc 1/x 2 (the other singularity at x = —xb is outside the 
kinematic region since if > 0). This is consistent with if oc l/l^ — > at large which 
ensures that the independently hadronizing u-quark is nearly on-shell. Thus 

C(x B ,x)~^ [ —M*) (B.7) 
a; J z 

Using this and the expression (B.5) for the MPD in the cross section (4.7) we get 
da(ir + N -» -r* L X) {ee d g 2 C F ) 2 f dlf + ( f dz 



dM\ ~ QV(1 - x B )N c J 2-Kli 



(ii) Comparison with standard factorization 

We should compare the cross section (B.8) with a calculation [3] where a(ir + + N — > 
u + ^f* L {q) + X) is expressed in terms of the hard subprocess cross section a{Ti + {k) +d(l<i) — > 
u(h) + 72 (<?)) convoluted with fd/ P (lf /p + )- Parametrizing the momenta as 

h = (lfJT,l±) 

q = (q + ,q~,-l±) (B.9) 
h = (#,(T,0 ± ) 

we have l± ~ %MP~, if = if ~ l± = ^f^i an ^ — - The subprocess amplitude 
is in the BB limit 

rpf +r . * \ 4iee d g 2 C F if f dz mirrt 

r(7r + d -> n + 7 L ) = T - \ tt- / — <M^) (B.10) 



QV2iVc V ^ J z 
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where the signs correspond to the helicity ±5 of the «-quark. The subprocess cross section 

t < m n^Y (B , U) 



Q 4 iV c 7 2irl+ 
gives the full cross section 

f rJl + /+ 
a( ^+Ar -^ u + 1 l + X) = J ^(vr+d - « + it) f d/p (^) (B.12) 

Using tfZ^ = ^<? + = d-Mj^ Q 2 /[ s ^ _ (l — x b)\ this expression agrees with (B.8) derived from 
the MPD. 

(in) fdu/p /fi/p as ^Af ->• 00 loztfj x M P~ ^ 

Finally, we consider the contraction of the <i-quark fields in the MPD, which corre- 
sponds to 7r + iV — > d + 72 + X. Similarly to (B.l) we have 



(0|C(y 3 )^(0)|0) = / §^o(i^)e-^f a 



(B.13) 

l-=l 2 /i+ 



After the integrals over d 2 l± and dy% and defining y = y 1 — y 2 — y 3 we get (with a minus 
sign from reordering the fermion fields) 

f d u/ P = -47^53 / dy^dy^dy" j ~ l + 9{l + ) (N(p)\My-) 1 + M0)\N{p)) (B.14) 

x exp [-Ifo- - y 2 ~)(Z+ + £+) - 1(2/7" + y^)(l+ - - ifo - % - y-)(q + + Z + ) 

The integrals over set = and l + = —if — q + = —/J, giving 

/du/p = S(lf ~ ^ 0(-Z 2 + )/ s(p (-^ /p+) (B.15) 
where the antiquark PDF is 

fu/ P (-lf/p + ) =-^J ^exp(-^y-Z+)(iV(p)|^(y-) 7 +^(0)|iV(p)) (B.16) 

Now the dominant contribution of C(xb,x) in the cross section comes from its pole at 
I2 = 0, giving ff oc e„ as expected. 
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